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1. Introduction 

Integer Codes are codes defined over finite rings of integers modulo to, m G 
Z. Since these codes have a low encoding and decoding complexity, they had a 
significant development over the last years and are suitable for application in 
communication systems (see [Ko, Mo, li, Ha, Ma; 10]). 

Some other codes similar to the Integer Codes, such as for example codes 
over Gaussian integers ([Hu; 94], [Gh, Fr; 10], [Ne, In, Fa, Pa; 01], [Ri; 95]) or 
codes over Eisenstein-Jacobi integers, have been intensively studied in recent 
years. 

In this paper, we will extend the study of Integer Codes to codes over subsets 
of real algebras obtained by the Cayley-Dickson process. This idea comes in a 
natural way, starting from same ideas developed by Hubner in [Hu; 94], in which 
he regarded a finite field as a residue field of the Gaussian integer ring modulo a 
Gaussian prime, ideas extended to Hurwitz integers in [Gu; 13] and to a subset 
of the Octonions integers in [FI; 15]. In this way, we regard a finite field as a 
residue field modulo a prime element from V, where V is a subset of an algebra 
At (K), where At (K) is a real algebra obtained by the Cayley-Dickson process 
and V has a commutative and associative ring structure. In this way, we obtain 
an algorithm, called Main Algorithm, which allows us to find codes with a good 
rate. This algorithm offers more flexibility than other methods known until now, 
similar to Lenstra’s algorithm on elliptic curves compared with p—1 Pollard’s 
algorithm. It is well known that for a prime p, Lenstra’s algorithm replace the 
group Z* with the group of the rational points of an elliptic curve Ci over Zp 
and, if this algorithm failed, the curve will be replaced with another curve C 2 
over Zp and we can retake the algorithm (see [Si, Ta; 92]). 
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In the case of Main Algorithm, the algebra At (M) and w offer this kind of 
flexibility since, for the same prime p, these can be changed and the algorithm 
can be retake. 


2. Preliminaries 


In the following, we shortly recall the Cayley-Dickson process for the real 
algebras. Let A be a finite dimensional unitary algebra over the real field R 
with a scalar involution : A ^ A, a ^ a, i. e. a linear map satisfying 
the following relations: ab = ba, a = a, and a + a, aa G M for all a,b G A. 
The element a is called the conjugate of the element a, the linear form t : A —>■ 

R, t (a) = a + a and the quadratic form n : A —>• R, n (a) = aa are called 
the trace and the norm of the element a. Since n (a) = aa = a (t (a) — a ), 
it results that — t (a) a + n (a) = 0, for all elements a G A, therefore A is a 
quadratic algebra. 

Let 7 G R be a fixed non-zero element. On the vector space A© A, we define 
the following algebra multiplication 

( 01 , 02 ) ( 61 , & 2 ) = (ai^i + 762^,0162 + ^ 102 ). 

We obtain an algebra structure over A © A, denoted by (A, 7 ) and called the 
algebra obtained from A by the Cayley-Dickson process. We have dim (A, 7 ) = 

2 dim A. 

Let X G (A, 7 ), X = ( 01 , 02 ). The map 

“ : (A, 7 ) -)> (A, 7 ) , X ^ X = (oi, - 02 ), 

is a scalar involution of the algebra (A, 7 ), extending the involution ~ of the 
algebra A. 

If we take A = R and we apply this process t times, t > I, we obtain an 
algebra over R, 

In this algebra, the set {cq = 1, ei,..., e„_i}, n = 2‘, generates a basis with the 
properties: 

ej = 7 A, 7 i G R, 7 i 0, i G {1, ...,n - 1} 

and 

^iCj — CjCi — /dijCk^ (dij G R, jdij 0, i 7 ^ j, i^j G { 1, ...n I}, 

Pij and Cfc being uniquely determined by Ci and Cj. 

Algebras A* (R), obtained by the Cayley-Dickson process, are power-associative 
(i. e. the subalgebra < a; > of At (R), generated by any element x G At (R), is 
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associative), flexible (i. e. x{yx) = {xy)x = xyx, for all x,y € At (R)) and in 
general it is nonassociative. 

For t = 2 and 71 = 72 = — 1, we obtain the Quaternion division algebra, 
Q (R), for t = 3 and 71 =72 = 73 = —1, we obtain the Octonion division 
algebra, O (R), and for t = 4 and 71 = 72 = 73 = 74 = — 1, we obtain the 
Sedenion algebra, S (R). Due to the Hurwitz’s Theorem, for t > 4, all obtained 
algebras are not division algebras (i.e. we can find a,b G At (R), a 0, 6 ^ 0, 
such that ab = 0 ). 

Let B = {1, 62 , 62 *} be the basis in Af (R), where 1 is the unity. If a; = 
2 * 2 * 

xi + J^XiCi € At (R), then its conjugate is the element x = xi — ^XiCi and 

i=2 i=2 

2* 

the norm of the element x is n (x) = xx = xx = '^xf. The norm n, in general, 

i—1 

is not multiplicative, i.e. for x,y G At (R), we have n (xy) ^ n (x) n (y). The 

2 * 

real part of the element x is Xi and its veetor part is '^XiCi G At (R). 

i=2 

2* 

The set At (Z) = {z G At (R) | z = xi + ^xtCt,xt G Z,i G {2,3,..., 2*}} 

i=2 

is called the integer elements of the real algebra At (R). This set has a ring 
structure, (see [Ma, Be, Ga; 09]) 

2 * 

Let w = a(l + X) eQ G At (R), a G R, and let V = {a + fere | a,b G Z}. We 

i=2 

note that t (x) = 2a, n (x) = 2‘a^ and — 2aw + 2‘a^ = 0. Since the algebra 
At (R) is a power associative algebra, it results that V is an associative and a 
commutative ring. 

Remark 2.1. For x G V, we know that the following properties are equiv¬ 
alent: 

i) X is an invertible element in the algebra V. 

ii) n (x) = 1. 

iii) X G {±1}. 

An element x G V is a prime element in V if x is not an invertible element 
in V and if x = ab, it results that a or 6 is an invertible element in V. 


Proposition 2.2. i) For x, y G V, we have n {xy) = n (x) n {y). 
ii) The ring Y is a division ring. 


Proof, i) Denoting with q = 2* — 1, let x = a + bw and y = c + dw. We 


obtain 


n (x) n {y) 


(a + ba)^ + b^a^q 


(c + da)^ + fl^a^q 


= {2aba a?' + b^a^ -I- b^qa^) {2cda + d?a^ -I- = 2ahc^a+2a^cda + 

Aabcda^ + + 2abd‘^a^ + 2b'^cda^ + 2abd^qa^ + 2h‘^cdqa^ + -b b^c^a"^ + 

+ a^fl^qa^ -b b'^c^qa^ -b 2b’^d‘^qa‘^ + b'^d^q^d^. 

Computing n {xy ), we get 
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n (xy) = [ac + {ad + be) a — a^bd (q + 1) + 2a^bd] ^+qa^ [ad + bc + 2abd\^ = 
2abc^a + 2a?cda + Aabeda^ + a^c^ + 2abd?a^ + 2b^cda^ + 2abd?qa^ + 2b'^cdqa^ + 
+ b'^c^a^ + b'^d'^a^ + a^d^qa^ + b'^c^qa^ + 2b‘^d'^qa‘^ + b'^d'^q^a^. 
Therefore n {xy) = n (x) n {y). 
ii) It results from i). 


Remark 2.3. The above result is also true for all elements from the set 
V'= {a + bw \ a,b € M}. 

In the following, we will consider a = -^,r>t — l,t>2. 

Proposition 2.4. If x,y £ Y, y ^ 0, with t > 2, then there are z,v gY 
such that X = zy + V and n (ri) < n {y). 


Proof. Since y 0 , we have that y is an invertible element in 
therefore | = a + bw, a, b 
\b — n\ < i. For z = m ' 


nw 


G R. Let m,n € Z such that |a — m| < | and 
G V and V = y [{a — m) + {b — n)w] , it results 


that - = z + - , therefore x = zy + v and v = x — zy. From here, we have that 
w G V! We obtain that 




"-1 


n (ri) = n (y) n {{a — m) + {b — n)w) = 
= n(y) [{a - m) + ^ {b - n)]"^ + 

= n (y) = (y) < n (y) .□ 


{b — nY 


< { %itll + (y) = 


Definition 2.6. With the above notations, let tt = x+yw be a prime integer 
in V and vi, V 2 be two elements in V. If there is u G V such that vi — V 2 = vir, 
then vi,V 2 are called congruent modulo tt and it is denoted by vi = V 2 mod tt. 


Proposition 2.7. The above relation is an equivalence relation on Y. The 
set of equivalence classes is denoted by Y^ and is called the residue classes of 
Y modulo TT. 

Proof. Denoting the elements from V,r in bold, if vi = V 2 mod tt and V 2 = 
mod TT, then there are v,v' G V such that vi — V2 = vtt and V2 — V3 = v'tt. It 
results that vi — V 3 = {v + v')tt, therefore the transitivity holds. □ 

Proposition 2.8. For each x,y G V, there is 6 = {x,y), the greatest 
common divisor of x and y. We also have that there are 7 and t G V, such 
that d = 7 a; + Ty.(the Bezout’s Theorem). 

Proof. We denote by J = { 7 X + ry | 7 ,t G V}.We remark that if 2 ; = 
'y'x + r'y G J and w £ Y, we have wz = {wj')^ + {wT')y £ J- We consider 
(5i = 71 a; + Tiy G J, such that (5i has the norm n (5i) minimum in J. We will 
prove that d = di.From Proposition 2.4, we have that x = + ri, with 

n(ri) < n((5i),yi,ri G V and ri = x — qi5i G J. Since n(ri) < n((5i) and 
5i £ J has minimum norm in J, it results that ri = 0, therefore di | x. In 
the same way, we will prove that (5i | y. Since di = 71 X + riy, it results that 
each common divisor for x and y is a divisor for 5i, therefore 5 \ di and finally 
(5 = 5i.n 


The above proposition generalized to elements in V Proposition 2.1.4. from 
[Da, Sa, Va;03]. 
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Proposition 2.9. is a field isomorphic to Jjjp'L^ p = n(7r),p a prime 
number. 

Proof. 

For Vi, V 2 € we define V 1 +V 2 = (fi + V 2 )mod tt and vi-V 2 = {viV 2 )mod 
TT. These multiplications are well defined. Indeed, if vi = v[ mod tt and V 2 = 
V 2 mod TT, it results that vi — v'l = un, V 2 — V 2 = u'tt, u, u' € V, therefore 
(vi + V 2 ) — {v'l + v' 2 ) = {u + u') TT. Since vi = v'l + utt, V 2 = V 2 + u'tt, it results 
that V 1 V 2 = v'iv '2 + M,r, with a multiple of tt. 

Denoting in bold the equivalence classes from Zp, let / be the map 

/ : Zp — )■ V^, / (m) = (m + tt) mod tt, where m S m. (2.1) 


Map / is well defined. Indeed, if m = to' mod p we have (to + tt) — (to' + tt) = 
m — m' = pq = TTirq, <7 € Z, therefore (to + tt) = (to' + tt) mod tt. 

From Proposition 2.8, we have 1 = Ui7r + U 2 ^. If / (m) = v,v = (to + tt) mod 
TT G Ytt, we define f~^ (v) = to (uitt) + to (v 2 ^) = to. 

Map / is a ring morphism. Indeed, / (ni) + / (m') = (to+ 7r) mod tt + 
(to' + TT)mod TT = (to + to' + TT)mod TT = f (m + m') and 
/ (m) / (m') = (to + tt) (to' + TT)modTT = 

= (toto' + (to + to') tt + Tr'^^mod tt = {mm' + TT)mod tt. We obtain that is 
isomorphic to Zp.D 

Let X = a + bw &Y, therefore we have n (x) = {a + ba)^ + q{ba)'^. For 
g = 2* — 1 and for certain values of t, we know the form of some prime numbers, 
as we can see in the proposition below. 

Proposition 2.9. ([Co; 89], [Sa; 14]) 

Let p gN be a prime number. 

1) There are integers a, b such that p = a'^ + 36^ if and only if p = l{mod 
3) or p = 3. 

2) There are integers a, b such that p = af+7b‘^ if and only if p = 1,2, A{mod 
7) or p = 7. 

3) There are integers a, b such that p = of + 156^ if and only if p = 
1,19,31, ‘i9{mod 60).□ 

The label Algorithm for At (R). 

1. We will fix t, a and therefore w. 

2. We consider tt G V a prime element, tt = a + bw,a,b G Z, such that 
n (tt) = p = (a + ba)^ + q{baY, with p a prime positive number. 

3. Let s G Z be the only solution to the equation a + &a; = 0 mod p, x G 
{0,l,2,...,p-l}._ 

4. Let r = G N, where [ ] denotes the integer part. 

5. Let /c G Z and k GZp be its equivalence class modulo p. 

6. For all integers a,T G { — r — l,...,r}, let c = (sr + cr) mod p and 
d = {a + raY + q{TaY. 
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6 . It d < p and c = k, then we find the pairs (cr, r) such that k is the label of 
the element a+rw € Vtt, that means a+rs = k mod p and n (ct + tw) is minimum. 
If we find more than two pairs satisfying the last condition, then we will choose 
that pair with the property that \<j\ + |t| < |a| + |6|. If there exist more than two 
pairs satisfying the last inequality, then we will choose one of them randomly. 


3. Codes over 


In the following, we will recall some definitions, which will be used in this 
section. 

We consider the ring of Gaussian integers, Z[i] = {a + hi | a, 6 € Z, = —I}. 
We know that a prime integer p of the form p = 1 mod 4 can be written of the 
form p = ttW, where Tr,W G Z[i] and W is the conjugate of tt. Let (Z[i]),n. 
be the set of the residue classes modulo tt. A block code C of length n over 
(Z[z])^ is defined to be a set of codewords c = (ci,...,c„), where Ci G (Z[z]) 7 r, 
i G {I, 2, ..., n}. For a', /3', 7' G (Z[z])^, with 7' = a' — /3' mod tt, the Mannheim 
weight of 7', denoted by wm (7O j is defined as 

WM (7O = \Re{'l')\ + , 

where Re{'y') represents the real part of the element 7 ' and Im{Y) represents 
the imaginary part of the element Y. Using the Mannheim weight, we can define 
the Mannheim distance between a' and j3' , denoted by du, as follows 

dM ia',l3') = WM ( 7 O • 


For other details, the readers are referred to [Hu; 94] . 

Using ideas from the above definition and generalizing the Hurwitz weight 
from [Gu; 13] and Cayley-Dickson weight for the octonions, from [FI; 15], in 
the same manner, we define the generalized Cayley-Dickson weight, for algebras 
obtained by the Gayley-Dickson process, denoted da- We will fix t, a, w, there¬ 
fore we will consider the elements in the algebra At (M). Let tt be a prime in V, 
TT =a -b bw. Let cr G V, cr = oo -b bow. The generalized Cayley-Dickson weight of 
X is defined as wg (x) = |ao| -b |6o| , where x = oq -b b^w mod tt, with jaoj -b j^oj 
minimum. 

The generalized Cayley-Dickson distance between x,y G is defined as 

da {x, y) = wg{x- y) 

and we will prove that da is a metric. Indeed, for x,y,z G Vjr, we have 
do (x,y) = WG (ai) = |ai| -b |&i|, where ai = x — y = ai -\- biw mod tt is 
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an element in V,r and |ai| + |&i| is minimum. 

da (y, z) = WG (<22) = |a2| + 1^21, where a2 = y — z = 02 + 62 v/mod tt is an 
element in and |a2| + I&2I is minimum. 

da (x, z) = WG (as) = losj + I63I , where = x — z = as + b^w mod tt is an 
element in Vtt and losj + I&3I is minimum. 

We obtain x — z = ai + a2 mod tt and it results that wg (oi + 02) > WG (<a3), 
since wg (03) = losj + I63I is minimum, therefore dc (x, y) + dG (y, z) > dG (x, z). 


In the following, we assume that tt is a prime in V with n (tt) = p a prime 
number of the form n (tt) = Mn + 1, M, n € Z, n > 0, such that there are /3 
a primitive element (of order p — 1) in V^n-, with the properties = w or 

= —w. We will consider codes of length n = 2^. 

The below definitions and Theorems adapted and generalized to all algebras 
obtained by the Cayley-Dickson process some definitions from [Gu; 13], [Ne, In, 
Fa, Pa; 01], [FI; 15] and Theorems 7,8,9,10,11,13,14,15 from [Ne, In, Fa, Pa; 01], 
Theorems 4,5,6,7 from [Gu; 13] and Theorems 2.3, 2.5, 2.7, 2.9 from [FI; 15]. 

Let G be a code defined by the parity-check matrix H, 


( 1 
1 

V 1 


p 

pM+1 

I^Mk+1 


p2{M+l) 


^2(Mfc-|-l) 


pn-1 X 

p{n-l){M+l) 

p{n-l){Mk+l) j 


(3.1) 


with k < n. We know that c is a codeword in C if and only if idc* = 0. If we con- 

n —1 

sider the associate code polynomial c{x) = ^ we have that = 

i=0 

0 ,1 € {0,1,..., k}. For the polynomial p (x) = {x — P) {x — P^^^) ... (x — , 

since the elements P,P^~^^, are distinct, from [10], Lemma 8.1.6, we 

obtain that c (x) is divisible by g (x), the generator polynomial of the code C. 
Since p (x) / ( x” ± w), it results that G is a principal ideal in the ring / 

(x" ± w). 

If we suppose that a codeword polynomial c (x) is sent over a channel and 
the error pattern e (x) occurs, it results that the received polynomial is r (x) = 
c (x) -I- e (x). The vector corresponding to the polynomial r (x) = c (x) -I- e (x) is 
r = c -|- e and the syndrome oi r is S = Hr*, where H is the above parity-check 
matrix. 

Theorem 3.1. Let C be the code defined on by the parity check matrix 

H={ 1 p ... ) . (3.2) 

It results that, the code C is able to correct all error patterns of the form e (x) = 
ejX*, with 0 < WG (e*) < 1 . 

Proof. Let r (x) = c (x) -I- e (x) be the received polynomial, with c(x) 
the codeword polynomial and e(x) = e^x* denoting the error polynomial with 
0 < WG {ci) < 1. Since /3” = w or /3” = —w, it results that e* = We have 
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the syndrome S = ^ with i, L € Z, 0 < t, L < n — 1. If we reduce L 

modulo n, we obtain i, the location of the error, and from here, I = and 
the value of the error. □ 

Theorem 3.2. Let C he a code defined by the parity-check matrix 
1 fi /3^ ... /3"-i 

2 ^ pM-ei ^2(M+l) ^(n-l)(M+l) 

Then C can correct error patterns of the form e{x) = eiX^, 0 < z < n — 1, 
with Ci S 

Proof. We consider the received polynomial, r {x) = c (x) + e (a;) with c{x) 
the codeword polynomial and e{x) = Cja;* the error polynomial with Cj G V.^.. 
It results that the corresponding vector of the polynomial r(x) is r = c + e. We 
will compute the syndrome S of r. We have Ci = fifiO < j < Mn — 1. Therefore 
the syndrome is 




S=Hr*= 


Si = pi+3 = \ 


We obtain with z + j = Mi mod{p — 1) and p(M+i)i+j-M 2 . _ 

with (M + l)z + j = M 2 mod{p — 1). We get Mi = {M 2 — Mi) mod{p — 1), if 
there is, then the solution of the system is z = mod n and j = [Mi — i) 

mod{p — 1). From here, we can find the location and the value of the error.□ 

Theorem 3.3. Let C he a code defined by the parity-check matrix 

/ 1 /3 /32 ... /3"-i \ 

H = \ 1 y 32 (M+l) ^(n-l)(M+l) I ^ 

\ ^ p2M+l ^2(2M+1) ^ _ ^(n-l)(2M+l) j 

Then C can find the location and can correct error patterns of the form e {x) = 
CiX*, 0 < z < rz — 1 , with e* G or can only correct error patterns of the above 
mentioned form. 

Proof. Using notations from the above Theorem, we have Ci = fifiO < j < 
Mn — 1. Therefore the syndrome is 


S=Hr*= 


Si = /?*+'? = 

SM+I = /3(^+i)*+^- = 
S2M+1 = /3(2M+l)i+J' = /3^3 


/ 1 /3 /32 \ 

Since the matrix 1 ^2(m+i) j^g determinant equal to 

I ^ y32M+l ^2(2M+1) j 

0 ipM (^ 02 M _ 2 )^ ^ 0, it results that the rank of the matrix (3.4) is 3, then this 
system always has a solution. We obtain = 1 , with i-\-q = Mi mod{p — 

1 ), _ 2 ^^ with (M + 1) z + j = M 2 mod{p — 1), ^(2M+i)i+2-M3 _ 



1, with (2M + 1) i + j = M 3 mod{p— 1). We can find the location of the error if 
Mi = (M 2 — Ml) mod{p — 1) and Mi = (M 3 — M 2 ) mod{p — 1) or, equivalently, 
* = mod n = mod n and the value of the error if 

(Ml — i) mod{p — 1) = (M 2 — (M + l)i) mod{p — 1) = (M 3 — (2M + l)i) 
mod(p — 1)(= j). □ 

Theorem 3.9. Let C be a eode defined by the parity-check matrix 


( I P 

I pM+l p2(M+l) 

I p2M+l p2(2M+l) 

^ I pd.M+1 p2(3M+l) 


\ 

p{n-l){M+l) 
^(«-l)(2M+l) 
^(n-l)(3M+l) ^ 


(3.5) 


Then C can correct error patterns of the form e (x) = Cix'^-i-CjxP 0 < i,j < n—1, 
with Ci , Cj € Vtt ■ 


Proof. We will prove this in the general case, when we have two errors. 
We have e* = /3®' 0 and Cj = P* 7 ^ 0, q', t' G Z. We obtain the syndrome: 


/ 


S=Hr*= 


Si = 0*+® 


yj+t' 




SM+1 = + a(“+i)j+‘' 

S2M+1 = Q;(2M+l)i+g' ^{2M+l)j+tr 

V S3M+1 = a(3M+l)*+9' ^ „(3M+l)i+i' J 


We denote = A and = B and we get 




S=Hr*= 




(3.6) 


s^= A-\- B 

sm+i=P^^A + P^iB 

S2M+1= P^^^A + P^^^B 

V S3M+1=P^^^A + P^^^B j 

If the system (3.6) admits only one solution, then code C can correct two 
errors. First, we will prove the following Lemma. 

Lemma. With the above notations, if we have two errors, we obtain P^'^ 7 ^ 
P^PQ <i, 3 <n-l and S 1 S 2 M +1 ^m+i- 

Proof. If = P^P then = 1 and Mn / M{i—j), which is false. 

Supposing that S 1 S 2 M +1 — s\i+i = 0, we have S 1 S 2 M +1 = Sm+i- For x = /?*+'? , 
we obtain that P'^^'‘six + P'^^^s\ — P"^^^six = + P'^^^s\ + 

2P^^ (/3^* — P^^) six. It results {P^^ — P^^f^x^ + 2/3'^*+^^sia; — /3^^*six — 
p2Mjg^x = 0. From here, a; = 0 or a; = ^ ^ = si. If 

have X = = si, this implies = 0, which is false. 


we 


We now return to the proof of the Theorem and we know that the following 
conditions are fulfilled: /3^* 


^ P^PO <i ,3 <n-l and S 1 S 2 M +1 ^ s\[+i- For 


B = si — A, it results that 
A [pMi _ ^ 
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^ (/ 32 M. _ p2Mj^ ^ 

A - /33Mj) = S3M+1 - Si/3^^A 

We obtain S2M+1 - = (sm+i - si/ 3 ^^) (/3“*+/3“^) and S3M+1 - 

If we denote by sm = /3^* + and pm = , we have 


S2M+1 ~ sm+iSm +PmSi — 0 


and 


(sm+ 1 — si/3 ^'^) (sm — Pm) = S3M+1 — si 0 ^^K 
It results that 


Pm = 


SM+ISM — S2M+1 
Si 


and 


Sm(siS2M+1 — = S1S3M+I — SM+1S2M+1- 

Therefore, we obtain 


Sm = 


S1S3M+1 — SM+1S2M+1 
2 

S1S2M+I — S^_|_]^ 


Pm = 


SM+1S3M+1 — S2M+1 


SlS2A^+l Sjy^_j_-j^ 

From here, by solving the equation — smx + pm = 0, we find the locations 
and the values of the errors. □ 


4. Main algorithm and some examples 

Definition 4.1. The rate of a block code is i? = where k is the dimension 
of the block code and n is the length of the codewords. 

The Main Algorithm 

Let p be a prime number. 

1. We find a, 6, t G N such that we can write p under the form 

p = + {2* -l)b^. (4.1.) 

We remark that the values for a,b,t, if there exist, are not unique. Let 
{ai, bi, ti}, I G {1, 2,..., u} all solutions of the equation (4.1). 

2. Let p = rijMj + 1, with nj,Mj not unique such that rijMj = p — l,j G 
{1,2, ...,u}. 

3. For I G (1, 2,..., u} and for j G {1, 2,..., 1 ;}, we find the algebra At, (R), the 

2*i 

element w = ^r)-i (1+ X)®*) ^ (R) ,n > ti — l,Y C Af, (R), the element tt G 

i=2 
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V, such that n (tt) = p, we find V,r such that Vtt is isomorphic to Ip and we 
find /? e such that = w or = —w. 

If the elements {aj,6/,t/} don’t exist, then the algorithm stops. 

If we have at least a solution for the equation (4.1) but we don’t find for 
j e {1,2, ...,v} the element /3 € such that = w or = —w, then the 
algorithm stops. If we have solutions in both cases, then we go to the Step 4. 

4. For each solution (o;, b}, I G (1, 2,..., u}, let J C {1, 2,..., v}. For each 

j G J, we have tij such that = w or /3”^' = —w. We can change w by 
increasing the value of r/ , if it is necessary, but working in the algebra At, (R). 
For each nj we compute Mj and the rate of the obtained code, Rj = Since 
we can suppose that the obtained codes have the same dimension k = kj, we 
will chose the indices I G (1,2,..., u}, j G J, the pair jo/, and the number 
Uj such that the rate Rj has the biggest value. 

In the following, we will denote by Algorithm 1, the method described in 
[Gu; 13] and by Algorithm 2, the method described in [FI; 15]. 

Remark 4.2. In the papers [Gu; 13] and [FI; 15] were developed some 
algorithms which built binary block codes over subsets of integers in the real 
quaternion division algebra and in the real octonion division algebra. The above 
algorithm generalized these two algorithms to real algebras obtained by the 
Cayley-Dickson process. Moreover, the Main Algorithm can be generalized to 
almost all prime number, which in general the Algorithm 1 and the Algorithm 
2 don’t make it. That means, in general, for a prime number p, we can get 
the algebra Aj (R), the element w £ At (R), the subset V C At, (R), tt G V, 
such that n (tt) = p, we can find V.^- with isomorphic to Zp, such that the 
obtained binary block code can have the highest rate. 

With the above algorithm, we have a much higher flexibility, similar to the 
Lenstra algorithm for elliptic curves compared with p — 1 Pollard algorithm. 
It is well known that for a prime p, Lenstra’s algorithm replace the group Z* 
with the group of the rational points of an elliptic curve Ci over Zp and if this 
algorithm failed, the curve will be replaced with another curve C 2 over Zp and 
we can retake the algorithm (see [Si, Ta; 92]). 

In the case of Main Algorithm, the algebra At (R) and w offer this kind of 
flexibility since, for the same prime p, these can be changed and the algorithm 
can be retake, with better chances of success. 

We will explain this in the following examples. 

Example 4.3. Let p = 29. We have a = l,b = 2 and t = 3, therefore 
p = 1 + 7-4 with unique decomposition. It results that we can use the real 

Octonion algebra. If we apply Algorithm 2, we have w = | ^1 + > 

TT = — 1 + 4:W,p = 29, n = 4, s = 22, /3 = 1 — w, = —w mod tt, therefore we 
can define codes. 

If we apply the Main Algorithm for w = ^ ^1 + ^ , we have tt = — 1 + 

8w,n = 4, s = 11 which is the label for the element w G Vtt. We remark that 


11 



we can’t find /3 € such that (3'*^ = w, as we can see from the MAPLE’s 
procedures below. 

for i from -15 to 14 do for j from -15 to 14 do 
c := (ll*j+i)mod 29; d := ((7/4)*j)~2+(i+(l/4)*j)"2; 
if d < 29 and c = 11 then print(i, j);fi;od;od; 

0 , 1 

4, -2 


A := 8~{-l} mod 29; for a to 29 do 
b := a~4 mod 29; if b = 11 then print(a);fi;od; 

11 


But, if we increase a. we still work on octonions and we take w 


32 


i + E< 


\ i=2 

with the label s = 24. We obtain /3 = — 1 — w with the label 4 such that 
/3^ = ru,therefore we can define codes. In this situation, both algorithms can be 
applied with success. 


Example 4.4. Letp = 71 = 64+7-1, with unique decomposition. Therefore 
a = 8,b = l,t = 3. Then we work on real Octonion algebra. If we apply the 

8 \ 

l+^ei),7 r = 7 + 2w,p = 71, n = 10, s = 

i=2 } 

82,13 = 2 — 2w, = w mod 7r.(see [FI; 15]) 

If we apply the Main Algorithm and first we take w = | ^1 + 

have TT = 7 + 4w,p = 71, n = 10,s = 16 which is the label for the element 
w € V^.We remark that we can’t find /3 € Ytt such that = w (even if we 
increase the value of r, as in Example 4.3), as we can see in the procedure below. 


Algorithm 2, we have w; = ^ ^ 


A := -7*4~{-l}mod71; for a from 1 to 71 do b := a~{10} mod 71; 
if b = 16 then print(a);fi;od: 

16 


Therefore, the Algorithm 2 is better than the Main Algorithm. 

Example 4.5. For p = 31 = 6 • 5 + 1, we have p = 4 + 3 • 9 = 16 + 15 • 1, 
therefore t € {4,16} and we can use the real Quaternion algebra or the real 
Sedenion algebra. If we apply the Main Algorithm for sedenions, we have w = 

I ^1 + ■ We get TT = 3 + 8w,p = 31, s = 19 and we remark that we 

can’t use it for the sedenions since we can’t find /3 G such that (3^ = w. 
Therefore, we will use the Main Algorithm only for Quaternion algebra, which 
can be applied in this case. 
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Example 4.6. Let p = 61. We have that p = 4- 3- 5 + l = 1 + 60 = 
1 + 15 • 4 = 49 + 3 • 4, therefore t € {4,16} and we can use the real Quaternion 
algebra or the real Sedenion algebra. 

If we take p under the form p = 61 = 7^ + 3 • 2^, we use the real Quaternion 

algebra. For w = ^ ^1 + ^ , we get tt = S+dru, the label for w is s = 14, n = 

10(p = 6-10 + 1) and we have (3 = —4 + w, = w, as we can see from in below 
procedures: 

A := -5*4~{-l}mod 61; for a to 61 do 
b := a~{10}mod 61; if b = 14 then print(a);fi;od; 

14 10 17 26 29 30 30 31 32 35 44 51 


for i from -31 to 30 do for j from -31 to 30 do 
c :=(14*j+i) mod 61 d := (3/4)*j"2+(i+(l/2)*j)"2; 
if d < 61 and c = 10 then print(i, j)fi;od;od; 

-4, 1 
1, 5 
5, -4 


In this case, the rate code is Rx = where k is the dimension of 

the code, since we can’t find (3 such that ( 3 ^ = w or f 3 ^^ = w, for Mj | p — 1, j G 

{1,2, ...,u}. 

If we consider p under the form p = 1 + 15 ■ 4, we use the real Sedenion 

/ 16 \ 

algebra, we get n = 4 and for w = | I 1 + ] , we have tt = —I + 16w, 

the label for u; is s = 42 and /3 = 2 + 2w. In this case, the rate of the code 
is i ?2 = f and it is greater than i?i. We remark that we can use both 

algebras to define codes, but in the second case, we have chance to obtain a 
better rate.(The dimension k is considered the same, in both situations). 


A :=16~{-l}mod 61; for a to 61 do b :=a"{4}mod61; 
if b = 42 then print(a);fi;od; 

42 25 30 31 36 


for i from -31 to 30 do for j from -31 to 30 do c :=42*j+i mod 61; 

d := (15/64)*j~2+(i+(l/8)*j)~2; if d < 61 and c = 25 then printCi, j) 
fi;do;do; 

- 6 , 8 
-5, -8 
-2, 5 
- 1 , -11 
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2 , 2 
3, -14 
6 , -1 


Example 4.7. Let p = 151 = 4 + 3 • 49 = 16 + 15 • 9 = 6 • 25 + 1. 

We have t € {2,4} and will use the real Quaternion algebra or real Sedenion 

algebra. For w = 1 ^1 + , we have tt = —3 + 14w,n = 25 and s = 140, 

the label for w. In this case, we can’t find an element /3, such that = w, 
= w, = w,etc, as we can see in the procedure below. 

A:=-3*14~{-l} mod 151; for a to 151 do b:=a~25 mod 151; 
if b = 140 then print(a);fi;od: 

140 


But, as we remarked, the number p can be written under the form p = 

16 + 15 ■ 9 = 25 • 6 + 1, then if we take t = 4, we can use the real Sedenion 

16 \ 

1 + ^ Ci ) . We obtain tt = 1 + 24w, n = 6 and 

i=2 J 

s = 44, the label for w. We can find P, such that P^ = w mod tt and /3 = 3 — 3w, 
with the label s = 22. 

A:=-24~{-l}mod 151; for a to 151 do 

b:=a~6 mod 151; if b = 44 then print(a);fi;do; 

44 22 51 100 122 129 

for i from -76 to 75 do for j from -76 to 75 do 
c := 44*j+i mod 151; d:= (15/64)*j"2+(i+(l/8)*j)~2; 
if d < 151 and c = 22 then print(i, j);fi;od;od; 


-9, 

11 

-4, 

-20 

-3, 

, 4 

3, 

-3 

4, 

21 

9, - 

-10 


algebra. We consider w = -^ 


Example 4.8. Let p = 149 = 25 + 31 • 4 = 121 + 7 • 4. In this situation, 
t G {3, 5} and we can use the real Octonion algebra or a real Cayley-Dickson 
algebra of dimension 32. 

We can’t use the Algorithm 2 for octonions, since we can’t obtain the element 
P and p is not under the form 7fc+l. For w = i I 1 + ^ j , we have tt = 9+8w. 

\ i=2 / 

We consider p = 1 + 4 • 37 and we can’t find an element P, even if we take 
p = 2k + 1 or 4A: + 1 or 37fc + 1. 
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A 

:= -9*8~{-l} mod 149; for a to 149 do 


b 

:= a~2 mod 149; 

if b = 92 then print(a) 

;fi;od; 



92 


A 

:= -9*8~{-l}mod 

149; for a to 149 do 


b 

:= a~4 mod 149; 

if b = 92 then print(a) 

;fi;od; 



92 


A 

:= -9*8~{-l}mod 

149; for a to 149 do 


b 

:= a~37 mod 149 

if b = 92 then print(a) 

;fi;od; 


92 


But we can choose another a. For example, for ic = i ^ e^J , we have 

TT = 9 + 16w, and si = 46, the label for w. If we consider p = 74n +1, n = 2, 
we get P = —2 + 4w, with label S 2 = 33. In this case, the rate of the code is 
^1 = For P = 37n + I, n = 4, the label oi P = Aw is S 3 = 35. In this 

case the rate of the code is R 2 = = |. We have i ?2 < Ri- Therefore the 

code in the first case is better, since it can have a greater rate as in the second 
case. For p = 2fc + I or 4fc + I, we can’t find /3. 

A := mod 149; for a to 149 do 

b := a~2 mod 149; if b = 46 then print(a);fi;od; 

46 33 116 


for i from -75 to 74 do for j from -75 to 74 do 
c := (46*j+i) mod 149); d := (31/64)*j~2+(i+(l/8)*j)~2; 
if d < 149 and c = 33 then print(i, j);fi;do;do; 

-4, 17 
-2, 4 
0, -9 
9, 7 
11 , -6 

A := -9*16~{-1} mod 149; for a to 149 do 
b := a~4 mod 149; if b = 46 then print(a);fi;do; 

46 35 50 99 114 


for i from -75 to 74 do for j from -75 to 74 do 
c :=(46*j+i) mod 149); d := (31/64)*j~2+(i+(l/8)*j)'2; 
if d < 149 and c = 35 then print(i, j);fi;do;do; 

- 11 , 1 

-2, 17 
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0, 4 
2, -9 


A := mod 149; for a to 149 do 

b := a~37 mod 149; if b = 46 then print(a);fi;od; 

46 


/ 32 \ 

If we work on a real algebra of dimension 32, let w = I 1 + ^ ] . We 

have TT = 3 + ?>2w, s = 107, the label for m, /3 = 4, with the label s = 4, n = 
4,p = 37 • 4 + 1, as we can see in the procedures below. 

A := -3*32~{-l} mod 149; for a to 149 do 
b := a~4 mod 149; if b = 107 then print(a);fi;od; 

107 4 27 122 145 


for i from -75 to 74 do for j from -75 to 74 do 
c := (107*j+i) mod 149; d := (31/256)*j~2+(i+(l/16)*j)"2; 
if d < 149 and c = 4 then print(i, j);fi;od;od; 

- 8 , 21 
-7, -18 
-4, 14 
-3, -25 
0, 7 
1, -32 
4, 0 
8, -7 
12, -14 


We can work on both algebras to obtain codes with good rate. 

Conclusions. Regarding a finite field as a residue field modulo a prime 
element from V, where V is a subset of a real algebra obtained by the Cayley- 
Dickson process with a commutative ring structure, in this paper, we obtain an 
algorithm, called Main Algorithm, which allows us to find codes with a good 
rate. This algorithm offers more flexibility than other methods known until now, 
similar to Lenstra’s algorithm on elliptic curves compared with p — 1 Pollard’s 
algorithm. 

As a further research, we intend to improve this algorithm. This thing can 
be done if first we can find answers at the following questions: 

i) When the equation (4.1) has solutions? 

ii) When we have more than one solution for the equation (4.1)? 

iii) If we have solution(s) for equation (4.1), when we can find the element 

/ 3 ? 
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